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>>■ ABSTRACT 



in 



There have been several model analyses of the near and mid IR flux from the 
circumstellar ring around HR4796A. In one set of models, the 10 and 18 fim IR flux 
has been attributed to the reprocessing of stellar radiation by /im-size particles. Since 
^ . these particles are being blown away, on a dynamical time-scale, by the radiation 

pressure of IIR4796A, they must be continually replenished by the collisional fragments 
I of larger particles. If the ring appearance persisted for the life span (8 x 10^ yr) of 

I HR4796A, a parent-particle reservoir with a total mass greater than SOOM^ would be 

I needed. In order to avoid being conspicuous at longer wavelengths, most of the mass 

' must be contained in parent particles larger than 20 — 40 cm. In other models, it has 

■ been suggested that the IR flux from the rings is emitted by sufficiently large particles 
that survive the radiative blow out by their host star. In a gas free ring, greater than 

2 ■ 3.2/xm-size particles would survive radiative blow out and a total of W~'^M^ would 

^ ! be adequate to account for the observed IR flux. But, in the vicinity of a young star, 

^ ' the possibility that the dust ring is embedded within a residual protostellar gas disk 

cannot be ruled out. In a gas-rich environment, larger sizes (> 100;um) are needed for 
^ ' the particles to survive the radiative blow out. The total dust mass required to account 

' _ 1 

■ for the IR flux is < 10~ M,^. The combined influence of gas and stellar radiation may 
also account for the observed sharp inner boundary and rapidly fading outer boundary 
of the ring. The pressure gradient induced by a small (10%) amplitude variation 
in the surface density distribution of a low-mass gaseous disk would be sufficient to 
modify the rotation speed of the gas. The resulting hydrodynamic drag on modest-size 
(> 100/im) particles would be adequate to compensate for the turbulent stirring, 
radiative drag and radiation pressure such that they remain gravitationally bound to 
the system. The required surface density variation of the gas may be induced by 1) the 
perturbation of a low-mass planet or the binary companion HR 4796B, 2) the photo 
evaporation of the disk, or 3) from the variations in the viscous angular momentum 
transport and mass diffusion rate in the disk. We show that the structure of the dust 
ring is preserved during and after the gas is being depleted such that similar rings may 
be common among early-type stars. 
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Subject headings: circumstellar matter, planetary systems, stars: formation 

1. Introduction 

The observational discoveries (Smith 8z Terrile 1984) of /? Pic type dusty rings provide 
tantalizing hints of residual planctesimal disks similar to the Kuiper belt may exist around other 
stars (Artymowicz 1997). Perhaps the most intriguing dust rings are those around HR 4796A 
(Schneider et al. 1999) and HD 141569 (Weinberger et al. 1999). The sharp surface density 
variations near the ring boundaries of HR 4796A are remarkably similar to that found around 
Saturn's F ring (Cuzzi et al. 1984) and Uranus' e ring (Elliot &; Nicholson 1984). Under the 
action of viscous stress, these differentially rotating disks have a general tendency to diffuse in the 
radial direction (Lyndcn-Bcll & Pringle 1974). The expansion of planetary rings are prevented 
by the tidal torques of nearby shepherding satellites (Goldreich & Trcmaine 1979, 1982). The 
structural resemblance between circumstellar and planetary rings naturally leads to the conjecture 
that the former, too, may be perturbed by nearby hypothetical planets (Jura et al. 1995). The 
reported azimuthal brightness asymmetry (Telesco et al. 2000) provides a supporting evidence for 
the scenario that the ring may be gravitationally perturbed by a planetary or stellar companion 
(Wyatt et al. 2000). 

Another motivation for invoking the embedded planet scenario is the perseverance of the ring 
itself. The NICMOS image shows that the surface brightness has a maximum at 70 AU with 
a sharp decline at the smaller radii and an inner cutoff near 50 AU. Outside this region, the 
surface brightness decreases more gradually with radius to less than half of the peak brightness 
at 90AU (Schneider et al. 1999). The observed ratio of flux from the ring to that from the star 
increases with wavelength at 1.1 and 2.6 /xm, indicating the mean particle size is larger than these 
wavelength. But, using the flux ratio at 10 and 18 urn and assuming a Mie emissivity, Telesco 
et al. (2000) inferred a "characteristic" diameter for the dust particles to ~ 2 — 3/im which is in 
contrast to the lower particle-size limit of lOfim inferred from the 60/um radiation by Jura {et al. 
1993). The observed spectral energy distribution (SED) with wavelength A = 10— 10'^ //m (Koerner 
et al. 1998) has been modeled. Using single-size particles (> 100/xm), Augereau et al. (1999) 
were able to flt the observed SED between 10 — 10^/xm but not the resolved images of the ring. 
Both the SED and the resolved images were better modeled with reprocessing by particles with 
a power-law coUisional particlc-sizc distribution (Heller 1970, Mathis et al. 1977) with adopted 
minimum (amin) and maximum (amax) sizes of 10/im and a few meters respectively. In these 
models, since the models are applied to SED in the wavelength range of 10 — lO^iim, the actual 
values chosen for amin and Omax do not significantly affect the actual match. 

With the stellar luminosity = 35^0, small particles, if they exist, would be blown away 
by the radiation pressure on the time-scale of a few hundred years (Backman & Paresce 1993). 
In the models adopted by Augereau et al. (1999), all particles with sizes larger than Omin would 
survive the radiative blow-out effect. But for the models proposed by Telesco et al. (2000), the 
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"characteristic-size" particles would not be able to remain. In a companion paper, Wyatt et al. 
(2000) suggested that these small particles are probably the collisional fragment of a population of 
larger particles. In this replenishment scenario, it is still essential for the parent particles to have 
an adequate reservoir of total mass to replenish the /xm-size particles for HR 4796A's life span 
= 8 X 10^ yr (Stauffer et al. 1995, Jayawardhana et al. 1998). In order to be consistent with 
the observed SED and ring structure, most of these parent particles must also be located within 
70 ± 20 AU and have a surface density comparable to or less than that inferred by Augereau et al. 
(1999). One or more nearby planets may provide both the supply and confinement mechanism for 
the parent particles. 

An important issue to be resolved is the required mass for such perturbers. In HR 4796A, 
such confining planets must either form or migrate to the vicinity of the rings (which is well 
beyond the orbit of Neptune) within r*. In conventional theories, the first stage of planetary 
formation is the emergence of planetesimals and solid cores (Safronov 1969; Wetherill 1980). 
Only after these cores have acquired a few Earth-masses does dynamical accretion of gas become 
possible (Pollack et al. 1996). But, beyond the orbit of Neptune, the formation of a few Earth 
masses core requires at least several lO''' yr (> r^) unless the mass of solid material is much larger 
than that inferred from the minimum mass solar nebula model (Lissauer 1993). Assuming an 
initial surface density 10-20 times that of the minimum mass solar nebula (which corresponds to 
150-300 M0 within the annulus between ~ 58 — 82 AU), Kenyon et al. (1999) found it possible to 
form planetesimals larger than 10^ km within 10^ yr. They further suggested that gravitational 
perturbation by these large objects may excite large velocity dispersion that may lead to disruptive 
collisions. However, if this assumed surface density of the dust particles is augmented for the solar 
composition, the inferred mass of the gaseous ring 0.05 — O.IMq within 70ibl2 AU) would make 
it gravitationally unstable. Cameron (1978) and Boss (1998) have suggested that gravitational 
instability in the disk may indeed lead to the formation of giant planets directly from the disk gas. 
But, some numerical simulations show that the growth of non axisymmetric perturbations excited 
by gravitational instability are more likely to induce efficient transfer of angular momentum in the 
disk rather than fragmentation (Nelson et al. 1998, Burkcrt & Bodenheimer 1996). The initial 
surface density required for this conjecture is also large compared with the mass inferred for the 
outer regions of typical protostcllar disks (Bcckwith & Sargent 1993a). 

In this paper, we examine the persistence of the dust ring structure around HR 4796A. Since 
the host star is only 8 x 10^ yrs old, we hypothesize the existence of some residual gas analogous to 
those found around several young stellar objects with similar ages (Beckwith et al. 1990, Beckwith 
Sz Sargent 1993b, Zuckerman et al. 1995, Hartmann et al. 1998). An upper limit 1-7 of gas 
has been inferred from the mm continuum radiation by Greaves et al. (2000). (A larger gas mass 
may exist if a substantial amount of heavy element is depleted from the gas phase and is stored 
in particles larger than a few mm.) We explore the interaction between the dust, radiation, and 
residual gas in the disk. In §2, we briefly recapitulate the physical interaction between particles, 
stellar radiation, and residual gas. We show that the hydrodynamic drag tends to promote the 
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effect of radiation pressure and increase the critical particle size for radiative blow out. In §3, we 
provide some constraints on the size and total mass of particles in the ring around HR 4796A. We 
show that variations in the surface distribution of the gas provide a favorable condition for modest 
(more than 100 //m) size particles to congregate into a dust ring. This ring structure is preserved 
during and after the depletion of the disk gas. We suggest that these particles provide a supply 
of the dominant contributors to the IR emission of the dust ring. Finally in §4, we discuss the 
implications of these results on planet formation and the ubiquity of similar rings around other 
early type stars. 

2. The Dynamics of Residual Dust Particles 

In this section, we examine the interaction between stellar radiation, residual disk gas, and 
dust particles of various sizes. 

2.1. The gaseous disk 

For illustrative purpose, wc consider a simple model in which a disk of dust particles is 
embedded in a residual gaseous disk. The only free parameters in this model are those which 
describe the gas surface density as a function of the disk radius R. All other gas properties 

are deduced in a self-consistent way. 

The qualitative results presented here depend on the gross global gas distribution in the 
disk and are insensitive to the detailed expression of S(i?). In principle, the S distribution is 
determined by the efficiency of angular momentum transfer (e.g. via viscosity from turbulence) in 

the disk and the tidal perturbation by planetary or stellar companions. The former mechanism 
generally leads S to attain a power law function of R in most regions of the disk (Lin & 
Bodenheimer 1981; Ruden &; Lin 1986). The latter process leads to a localized S variation such as 
gaps (Lin & Papaloizou 1979; 1993) and confined rings (Goldreich & Tremaine 1978). Although 
low-mass companions may not be able to induce the formation of a gap, they can nevertheless 
introduce variations in the surface density of the disk. Self-gravitational effects will not affect the 
disk structure as the surface density in our models is always orders of magnitude smaller than any 
critical value (see Toomre 1964). For computational purpose (e.g. to produce our plots), we adopt 
the following generic prescription in which the global variation of S is non monotonic: 

if : R<Ri ^ 

if : Ri<R<Ro 
if: R>Ro J 

where d is a model parameter which determines the amplitude of a local maximum. The magnitude 
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of E at some fiducial radial location Rq is set to be Sq. For R < Ri and R > Rq, decreases 
monotonically with R following a power law (with an index -2.5) which is somewhat steeper than 
that for the minimum mass solar nebula model (Hayashi et al. 1985). But in the outer region of 
the disk, E is likely to decline more rapidly with R as inferred from observation (Mundy et al. 
1996, Hartmann et al. 1998). Since our results are insensitive to the detailed nature of the global 
radial dependence of S, we choose such a power law so that the disk mass remains finite at large 



The inner and outer regions of the disk are connected through a transition zone at 
Ri < R < Ro- In this zone the surface density has the shape of a Gaussian with a width of 
SRo. The value of Ri is set to satisfy {Ri - Ro) = 5Ro{2.5ln{Ri/Ro) - Ind)^^ such that S 
is continuous. In all our models, we choose d < 1 which is reasonable because the continuum 
flux in the mid-IR wavelength range is observed to be more intense than that in the near-IR 
range in several protostellar disks (Hillenbrand et al. 1998). One possible interpretation is that 
these disks are more depleted in the inner region or better preserved in the outer region as in 
the context of circumbinary disks (Jensen &; Mathieu 1997). On the theoretical side, such a E 
variation is expected if a) the viscosity decreases rapidly with radius (Ruden & Pollack 1991), b) 
the inner region is depleted by either a stellar or a disk wind (Shu et al. 1994, Konigl & Ruden 
1993), c) HR4796A is a AOV star and a strong source of UV radiation which could lead to photo 
evaporation in some regions of the disk (Shu et al. 1993) and quench gas infall along the rotation 
axis and onto the inner regions of the disk (Yorke et al. 1995), d) the disk is depleted by the 
formation of one or more planets which are expected to first emerge in the inner regions of the 
disk (Jura et al. 1993), or e) the structure of disk around HR4796A may be tidally perturbed by 
the stellar companion HR4796B (Augereau et al. 1999), especially if the binary orbit is eccentric 
(Korycansky &: Papaloizou 1995). The difference between the inner and outer regions of the disk 
can be minimized by setting d close to unity. Finally, the S variation in the transition region 
Ri < R < Ro can also represent a tidal perturbation of a nearby embedded planetary companion. 

The value of Sq is set such that the mass Mg = 27r/^^°+/^J' RJ:{R)dR contained within a ring 
centered on Ro with a half width SRo is to be specified as a model parameter. This mass is to be 
spread over a ring with a half width SRo = 10 — 20 AU centered on a radial location Ro = 70AU. 
In our models Mg is 1-100 Earth masses (Mq). The corresponding S is orders of magnitude 
smaller than that for the minimum mass solar nebula model at a similar distance from the central 
star (Hayashi, Nakazawa & Nakagawa 1985). 

The outer regions of the disk are optically thin such that every dust particle in the disk 
receives direct irradiation from the central object. The particles are heated to the local black 
body temperature (small grains may be hotter because they cannot radiate efficiently, Chiang 
& Goldreich 1997). The luminosity and mass of HR 4796A are = 35Lq and = 2.5M0 
respectively (Jura et al. 1993). The temperature of the dust grains can be given by: 



R. 




(2) 
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where a is the radiation density constant. This upper hmit estimation holds for particles larger 
than the typical wave length of the radiation and assuming an albedo of zero. A more careful 
treatment of the optical properties of the particles is not necessary in the context of the effects 
studied in this paper. 

If gas molecules collide frequently with the dust particles, they can also attain a similar 
temperature, T^igh = T^. In the low gas density limit, however, thermal conduction between 
the dust and the gas is less efHcient than the gas radiative processes. As a consequence of 
this thermal imbalance, the gas temperature may decrease to a minimum value {Tioy, = lOK) 
which is comparable to that of the radiation background field (Preibisch, Sonnhalter and Yorke 
1995; Yorke & Lin in preparation). In general, Tio^ < Tgas < Thigh and the associated sound 



speed Cs = yTgas * Rgas/fJ-^ where Rgas is the gas constant and the mean molecular weight 
(i.e. n = 2.353 for an He — H2 mixture). The disk thickness is determined by its pressure 
scale height Hp = Cg/O. where Q. = {GM^/R^)^/'^ is the Keplerian frequency. At i? = 70 AU, 
0.03 < Hp/R < 0.1 for the low and high temperature limits. For the parameters we have adopted, 
the disk gas is likely to be in the low temperature and small thickness range. The local gas density 
in the midplane of the disk can be estimated as p = H/lHp and the pressure p becomes p = c^p. 

The azimuthal orbital frequency of the gas (a;) is given by the condition for hydrostatic 
eqTiilibrium between the gravitational potential of the central object and the sum of centrifugal 
acceleration and the radial pressure gradient 

GM^ 9 1 

_l=^2^__Vp, (3) 

which leads, in first order approximation, to a local deviation dV from the Keplerian profile of 
(see Whipple 1972) 

,v.(.-a)R = v,-aR = ^^^,^A.^ (4) 

where = uR is the azimuthal speed. In §2.4, we consider the case that this angular frequency 
may be modified by the particles' drag when they become the dominant constituent of the disk. 

Results under specific assumptions for the free parameters of this model will be discussed in 
the §3. 



2.2. Radial drift of dust grains 



The particles are drifting radially under two independent physical effects: 1) hydrodynamic 
drag by the gas and 2) the radiation pressure of photons emitted by the central object. 
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2.2.1. Departure from Keplerian flow 

While the pressure gradient induces the gas to attain a finite dV, the particles' motion is 
determined by the gravity of the central star, the hydrodynamic drag by the disk gas, and the 
drag and pressure by the stellar radiation. Radial inward and outward drift of the particles can 
be generated by the non-Keplerian motion of the gas. Small particles are well coupled to the 
gas, thus they also move on non-Keplerian orbits provided their Stokes number St = TfQ. « 1.0 
(Weidenschilling 1977). The friction time for a spherical particle in the Epstein regime (e.g. Klahr 
& Henning 1997) is given by 

^ Csp 4:AdCsp' 

where ad is the particle's radius, its mass, = vra^ its cross section, pd its density, Cs the 
thermal velocity of the gas and p the local gas density. In the case of transonic velocities (which is 
irrelevant in our model) Cg has to be replaced by y/cj + dv'^ with the relative velocity (dv) between 
the particle and the gas (H. Yorke 1999, private communication). The critical particle size which 
demarcates the small and large particle size range then becomes 

ao = S/2pd (6) 

which is a few cm for a minimum mass solar nebula model at 70 AU and smaller for a depleted 
gaseous disk. In our standard models (see Sect. 3.5.1 and Sect. 3.5.2| ), ao = 0.06 cm at Rq. 



In the solar system today, the absence of gas makes St >> 1 for all particles. Modest-size 
particles spiral inwards due to the Poynting Robertson effect (see below), but spherical particles 
(with pd = 2.5gm cm~^) with size 

ad<ai = -^yrrr 'i-'ipm (7) 

iDvr GJvUcpd 

would be blown out by the radiation pressure on hyperbolic orbits as "beta meteorites" (see Griin 
et al. 1985) on a time-scale comparable to the local dynamical time-scale ~ il^^ ~ 10^ yr (see 
below). For our considerations we neglect particles smaller than the typical wavelength of the 
radiation (No Mie scattering). 

Larger particles would remain bound to the host star. They can even remain on a circular 
orbit with a sub Keplerian azimuthal velocity 



R 



(8) 



where denotes the particles theoretical velocity in the absence of gas, which is an important 
quantity (see below). While the particles are orbiting around the central star they also experience 
a radiative head wind due to the aberration of its photons. In gas-free environments, such as our 
own solar system, this Poynting-Robertson radiative drag effect induces the azimuthal velocity 



- 8 - 



{V(f)) of small particles to change at a rate 

'''^ = ~^^'^- 

This process leads an inward drift for all remaining particles with > ai. 

In a gaseous environment, the combined influence of hydrodynamic and radiative drag induces 
the particles' azimuthal velocity to change at a rate 

VS = ■ „n 7^V^. (10) 

If Tf is small in comparison to the relevant evolutionary time-scales (see Fig. || and Fig. an 
equilibrium would be rapidly established with = so that 

Substituting tj, we find 

which is independent of particle sizes (as long as 5t << 1 and the particles are larger than the 
wavelength of the radiation)! 

Around young stellar objects where the particles are embedded in a gaseous circumstellar 
disk environment, the term dv = v,f) — (see Fig. ||) is typically less than of 1 mm/s at 70 AU 
while dV ~ lO^cm/s. Note that the Poynting Robertson effect is insignificant until S is depleted 
below 3Lj, /Svr i2^c^r2 which is not only more than four orders of magnitude below that for our 
standard model and that inferred from mm continuum radiation for typical protostellar disks but 
also smaller than that inferred for the CO gas around some young stellar objects (Zuckerman 
et al. 1995). (The actual gas content may be even larger if a significant amount of CO gas is 
depleted due to grain condensation.) Thus, in our estimate of the particles' motion, we only 
consider contributions from dV and neglect the Poynting Robertson effect. 



2.2.2. Radial drift 

In addition to the radiative drag in the azimuthal direction, the stellar photons also exert on 
the particles a radiative pressure in the radial direction. In general, the radial drift Vr of small 
particles (with St < 1) results from the non-Keplerian motion and the radial radiation pressure 
can be obtained from the radial component of the equation of motion such that 

dvr _ Vr v<t>^ GNh Ad 

dt ~ Tf R i?2 + 4vri22^^c' ^ ' 
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In order to simplify work we define r2*(= v'^/R) which is the orbital frequency with respect to the 
radial radiation pressure (see Eq. ^). Note that fi* < in general. Thus it follows 

+ (14) 

dt Tf R 

Here dV* = V(f, — n*R denotes the deviation of the gas from the modified Keplerian orbit. 
Since the radius dependence of gravity and radiation pressure is identical, the outward radiation 
pressure effectively reduces the gravitational attraction on the particles imposed by the host star 
everywhere (Gustafson 1994). 

The "characteristic size" of particles inferred by Telesco et al. (2000) is smaller than ai 
such that they are expected to be blown away by the radiation pressure of HR4796A. In order 
to remove this dilemma, Wyatt et al. (1999) suggested that although these small particles are 
continually being blown out, they are also being replenished by collisional fragmentation of larger 
particles which are not significantly affected by the radiation pressure. In a gas free environment, 
the parent body size must be larger than ai. 

But, around many young stellar objects, small particles are embedded in a residual gaseous 
background. The gas generally does not directly experience the radiation pressure because its 
molecular opacity is too small. Thus, the gas azimuthal speed {V^) is independent of and it 
is given by eq. (^). The discussions in §2.2.1 indicate that particles with St < 1 (or equivalently 



with Qd < oo) have ~ (see eq. 12). As the hydrodynamic drag forces these small particles to 
nearly corotate with the gas, they attain a finite Vr as a consequence of being out of hydrostatic 
equilibrium. When a steady state is established, 

Vr = Tf2n*dV* for St«l, (15) 

in the limit that the quadratic term dv'^ is small in comparison to Q*Rdv* (cf Weidenschilling 1977; 
He used another definition of dVl). The sign of dV* can now be negative or positive, depending 
on whether the particle is less or even more sub-Keplerian than the gas. Relatively large particles 
follow in a similar fashion, 

Vr = dV* for St = 1, (16) 

and 

2dV* 

Vr = for St » 1. (17) 

Note that the effect of radiation pressure decreases with the particle size. From this 
equation one can deduce that in regions where dV > (as a positive of a negative gas pressure 
gradient), all the particles would migrate outwards since their azimuthal speed is always slightly 
hyper-Keplerian, just like the gas. But in the more common situation where dV < (which 
corresponds to a negative gas pressure gradient as expected in unperturbed regions of protostellar 
disks), hydrodynamic drag is more important than radiative drag for all particles with 

"'^ > " - 327rR^cnpddV ^^^^ 



- 10 - 



so that they undergo orbital decay. In contrast, particles with ad < 02 are blown away by the 
stellar radiation pressure (see Fig. Q and Fig. |2|) . This minimal size for the given parameters of HR 
4796A is 02 ~ 100 — 200/im for almost all disk radii R. Only at those location where variations in 
S lead to a reduction in dV is 02 greater than 500^m. In the location where dV vanishes or attains 
a negative value, 02 becomes infinite and particles of all sizes (including planets) are prevented 
from inward radial drift. 

We note that the magnitude of 02 is larger than oi in eq(0) by a factor 
~ {E? / Hp){d[np / d\nR)^^ 10 — 100 because the small particles are forced to corotate 
with the gas such that a centrifugal balance is always maintained. Thus, if they are embedded 
in a residual gaseous disk, the "characteristic-size" (^m-size) particles in HR 4796A inferred by 
Telesco et al. (2000) must be supplied by parent bodies which are at least two order of magnitude 
larger than them. 

In a gaseous disk environment, all particles with smaller than both qq and 02 are blown 
away by the radiation pressure on the dynamical time-scale r^. In the limit that oq > 02, particles 
with sizes oq > > 02 can achieve both positive and negative radial drift velocities depending on 
the local gas pressure gradient. In §3, we present numerical results to show that these changes 
in the direction of particles' radial drift lead to the particles' accumulation into rings. Note 
that the value of 02 does not directly depend on S. As the disk gas is depleted, neither dV nor 
a2 are strongly affected. But, oq is linearly proportional to S. When ao decreases below ad of 
some particles, they decouple from the gas and their azimuthal velocity becomes closer to the 
local modified Keplerian value (see eq. ^). Since their ad is larger than 02 initially and 02 
decreases, the radiation pressure and the Poynting Robertson effect have a weaker effect than the 
hydrodynamic drag throughout the evolution such that their radial drift velocity is described by 
eq. 17. Eventually, the depletion of the residual gas increases tj and reduces Vr- 



2.3. Constrains on the disk profile 

We already noted that the trapping mechanism for particles is located at places, where the 
disk becomes Keplerian, as the radial pressure gradient vanishes, i.e. 

Vp = V (c^p) = V {csm) = 0. (19) 

Since ~ the general criterion for the ring formation is 

cSoci?^/^ (20) 

This criterion can now be used to determine the local shape of the Surface density profile. We 
consider three distinct cases: 1) an isothermal disk(T = Tiow)-, 2) an optically thin disk in radiation 
equilibrium (T = Thigh)-, and 3) an optically thick viscous disk. 

In the isothermal case Cg is constant such that the trapping mechanism requires S to increase 
locally as i?^/^ regardless its magnitude and radial extent. In the standard optical thin case, 
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(perhaps most relevant to the case of HR4796A) , Cg ^ R 4 as T^i? 2. The necessary power-law 

7 

i?— dependence in S is only slightly steeper than in the previous case: i.e. S ^ i?4. 

In an optically thick disk which is heated by the local viscous dissipation, the disk's mid plane 
temperature = (3x^/8)^^^ where the vertical optical thickness ~ and the Rosseland 

mean opacity k = KqT^ (Lin & Papaloizou 1985). If the disk is approximately in a steady state, 
its effective temperature Tg cx (e.g. Lynden-Bell &: Pringle 1974, Bell et al. 1997) such that 

^ _R-3/25.i/2_ Particles are trapped in regions where T, on Kq ^^^R^/^ which is 
similar to the isothermal and the optically thin cases. 

The dominant contributors to the disk opacity are grains. Across some critical locations (dust 
destruction zone) in the disk where the grains sublimate, the opacity increases mainly through 
changes in the scaling constant Ko (Ruden &; Lin 1986). Particles may be trapped in regions where 
Ko oc R^. A full 2-D detailed model is needed to address the issue whether such a transition may 
indeed be a viable effect to prevent dust in form of meter sized solid bodies from drifting into the 
sun and help collect bodies at certain transition radii from the sun to form the terrestrial planets. 
For the present context, the ring around HR 4796 A is well outside the ice and silicate destruction 
zones. 

2.4. Particle-dominant flow 

In the above discussion, we neglect the feedback effect of hydrodynamic drag on the gas 
flow. In principle, the gas flow may be perturbed by the motion of the particles when the 

particles dominate in mass surface density (^d) over the surrounding gas (Nakagawa, Sekiya, & 
Hayashi 1986). In the Kuiper Belt today, no sign of residual gas has been detected. Gas may be 
preferentially depleted prior to the depletion of the dust particles. When the surface density of the 
gas (S) decreases below that of the dust (S,^), the hydrodynamic drag between them may induce 
the gas to comove with the particles rather than forcing the particles to corotate with the gas. In 
reality, however, the only particles which can provide an efficient momentum transfer with the gas 
are those with sizes « (Nakagawaet al. 1986). During the depletion of gas, oq decreases 
with S, leading to a reduction in the total fraction of particles which are well coupled with the 
gas. Thus, the gas speed is unlikely to be significantly affected by the particles' drag effect. 

For those regions of the disk with S > Ti^, gas flow would be affected by the hydrodynamic 
drag of the particles near the midplanc of the disk if the total mass density of the particle 
ndTTid exceeds the gas density p ~ ^/Hp (Weidenschilling &; Cuzzi 1993). (Similarly, particle 
accumulation due to variations in the radial drift speed may also lead to local enhancement of 
particle density.) The average mass density of the particles pdust is determined by their scale 
height Hd = CFd/^ where cj^ is the velocity dispersion of the particles. In a turbulent medium, the 
particles arc being stirred by the dispersive motion of the gas. Since the turbulent speed of the gas 
is unlikely to exceed its sound speed, particles with > ao would sediment toward the mid plane 
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and attain a scale height < H. With a substantiahy sub sonic turbulent speed, sedimentation 
is also possible for sufficiently small (a^ < ag) particles which couple well with the gas. Only after 
Pdust of the small (a^ < «o) particles has exceeded the gas density ~ T,/2H, would the gas motion 
be forced to corotate with the small particles' Q* of the rather than uj. Without a significant 
differential motion between them and the gas, particles with ai < Ud < oo would not undergo any 
further orbital migration. But larger (a^ > ao) particles are less affected by the radiation pressure 
and have a larger v^. Gas drag induced by the residual headwind would continue to drive the 
larger particles inward though at a much reduced pace. 

But, the optical depth of HR4796 A is well below unity and 

(see eq. below). For particles of = 200/im the density is Pdust « 6.3 x 10^^'^g cm~^, which is 
well below the gas density inferred from the minimum mass solar nebula model. For those models 
we have adopted in this paper, particles have negligible drag effect on the motion of the gas. 



3. Models for HR 4796A 

In this section, we present a scenario for the ring structure around HR 4796A. First, we revisit 
the two existing models for particle size distribution by analyzing the preferred size of particles 
for retention. We describe a simple quasi static model in which we assume that HR 4796A had 
initially a circumstellar gas and dust disk around it. We then show that a dust disk will evolve 
into a ring structure with sharp edges. We argue such a structure would be preserved after the gas 
is depleted. 



3.1. Replenishment of a population of pm-size particles 

First, we examine the implication of the scenario proposed by Telesco et al. (2000) in which 
the dominant contributors to the 10 and 18 nm radiation from the ring around HR 4796A are 
/im-size particles. Since the "characteristic" size of the particle is comparable to the wavelength 
of the photon, the reprocessed radiation by the disk particles 

Ldisk = QLi,- — ^ — (22) 
47r/t^ nid 

where Q is the albedo and assumed to be one for large particles and dR the radial width of the 
dust ring. Assuming SiO-dust particles which have spherical shape and a density pd = 2.5gm 
cm~^, the total mass of those particles with "characteristic" sizes is 

16tT Ldisk r>2 /r,o^ 
Mchar = ^ p R Pdttd- (23) 
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The total IR luminosity of the ring is estimated to be Ldisk/Li, = 5 x 10 ^ (Jura et al. 1993) so 
that 

Mchar ~ (fld/lcm) X 2.3 X lO'^^g ^ (arf/lcm) x 38M^. (24) 

For the all particle size to be at the "characteristic" size Oc = 1 — 2fim, one would need 
- 4-8 X lO-^Me. 

If the IR emitting particles are smaller than both ai and 02 they would be blown away from the 
ring region on a dynamical time-scale = ~ 10^ yr in both gas-rich and gas-free environment. 
Over the life span of HR4796A (r^ = 8 x 10^ yr), an amount M[^J « Mchar^r^ ~ 300Me of 
characteristic-size particles is lost through blow out by radiation pressure. In order to replenish 
this loss, Wyatt et al. (2000) proposed that the IR-emitting /xm-size particles are collisional 
fragments of larger parent particles. Unless we are observing HR 4796A in a special epoch, the 
total amount of mass in the parent particles (Mfot) is likely to be comparable to or larger than 

In the model presented by Telesco et al. (2000), only a single population of characteristic 
//m-size particles is assumed to be responsible for the reprocessing of stellar radiation. In order 
for the parent particles assumed by Wyatt et al. (2000) to remain inconspicuous, their total 
surface area must be less than that of the /im-size IR-emitting particles. For a population 
of parent particles to contribute a major fraction of Mtot, their total area is oc a^^ . Since 
^ tost' / char ~ ~ 10^, the minimum size for the parent particles is amm ~ lO^Cc ~ 10cm such 
that their Tz is less than that of the IR-emitting /xm-size particles. 



3.2. Particle collisions 

The production of /xm-size particles requires collisions. The collisional frequency for individual 
characteristic particles is 

Uc ~ UdAdad, (25) 

where dd is the velocity dispersion of the particles. The spatial number density of those particles 
which contribute to most of the IR radiation is 

_ Mchar ^ 4:R Ldisk Pdad 

md4.TrRdRHd~ 3dR rudHd 
For photons with wavelength smaller than the characteristic size of the particles, the particles' 
opacity is determined by their geometric cross section Ad- The optical depth in the direction 



normal to the plane of the disk is ~ {Mchar /fT^d) {Ad 1271 RdR) ~ Ldisk/ Li,. From Eqs. (23) and 
(p5| ) we find the collisional frequency for IR emitting particles to be 

~ = TB^^ ~ « ^d'- (27) 

art -L^ 

If the dominant contributors of the IR radiation are /xm-size particles as postulated by Wyatt et 
al. (2000), they would be blown out on a dynamical time-scale because their < ai (in a gas- free 
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environment) and aa < 02 (when the particles are embedded in a gaseous disk). Thus, after they 
are produced from the colhsions of their parent particles, these particles would not have sufficient 
time to collide, fragment, or coagulate among themselves before they are ejected of the ring region. 
Furthermore, in order for the ^m-size particles to remain as the dominant IR emitters, both 
and Vc must decrease with particle size so that the particles with 1 — 2ij,m < < ai or 02 would 
not dominate the radiation reprocessing. Since both ai and 02 are considerably larger than a 
few ^m, the hypothetical fim-size IR-emitting particles must be produced directly (rather than 
through a cascade process) during collisions between the smallest surviving particles (those with 
aa > ai or > 02). 

We now consider the collision frequency of the population of parent particles which are the 
main suppliers of the //m-size IR-emitting particles. In the minimum-size (with ~ Qmm ~ 10 
cm) limit, the parent particles' Vc is comparable to that of the IR-emitting particles, i.e. Vci in 
eq(p7|). On the dynamical time-scale during which the ^m-size particles are blown away, a total 
fraction I'ciTd of the parent particle population would collide. In order for this limited number 
of events to directly replenish the mass loss due to the radiative blow-away of the IR-emitting 
particles, each collision must generate Nf ~ {r^/Tiif'Li./Liiisk ~ 10^^ (/im-size) fragments. The 
parent bodies could have a size > amin, in which case, collisions would occur less frequently 
by a factor ~ amin/dd provided they can lead to the production of ~ {Ti,/Td)'^{adL^,/aminLdisk) 
fragments. Note that particles with sizes > {Ti,L(iisk/TdLi,)amin collide less than once during r^. 
Nevertheless, at least in principle, collisions among a fraction of the large-particle population may 
replenish the //m-size IR-emitting particles. 

The considerations presented in this subsection indicate that the /im-size particle scenario 
suggested by Telesco et al. (2000) and Wyatt et al. (2000) requires not only a population of parent 
bodies which are larger than a few cm but also a non continuous bimodal (~ l^m and > 10cm) 
particle-size distribution. There is no simple physical mechanism to preferentially produce this 
type of particle distribution function. 

3.3. Particle size distribution 

In reality, collisions are more likely to produce fragments with a range of sizes. As noted 
above that all fragments with > oi (in a gas-free environment) or > 02 (if the fragments are 
embedded in a residual gaseous disk) would survive the radiative blow out. During the life span 
of IIR4796A, the minimum-size parent particles would undergo VciTi, ~ {Ldisk / Li,){Ti, / Td) ~ 10^ 
collisions. Although fragmentary particles with oi < < amin or 02 < < amin may collide less 
frequently, they may nevertheless attain a collisional coagulation-fragmentation equilibrium. 

In principle, the particles' size distribution may be obtained by solving the appropriate 
coagulation equation with both coagulation and fragmentation effects included. Approximate 
solutions may be obtained for collisional steady state in which the disruption or growth rates {Md) 
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of the total mass (M^ = Ndrud) for a given population may be assumed to be independent of ad- 
If we assume the size-distribution evolution occurs primarily through collisions between similar 
size particles, 

Md ^ MdUc ~ Nimd{ad/Rfn. (28) 

A constant Md would imply Nd oc a^^'^, Md oc a^'^ and oc a^*^'^, i-e. most of the mass are 
contained in the large particles whereas most of the surface area cross sections are contained in the 
small particles. This particle size distribution is similar to that derived for collisional equilibrium 
(Hellyer 1970, Mathis et al. 1977) and adopted by Augereau et al. (1999) in which 

dNd/dad oc af-^. (29) 

More generally, the growth and fragmentation may be mainly regulated by collisions between 
different-size particles, which would result in a slightly modified size distribution. 

The distribution function in eq(|29|) is only applicable for particles with ^ ai or ~ a2. 
Since smaller particles are subject to radiative blow out of the ring before they can collide with 
each other again, their size distribution is probably determined by that resulting from the break 
up of the parent bodies. But if this distribution function can be extended to the small particles the 
total mass carried by the blown-away dust particles is Ml°l^ ~ Mfj^.f (oi/oc)^/^ ~ 5 - 10 x lO^M® 
in the gas free environment and M/°*^ ~ M/^'^^'"(a2/ac)^/^ ~ 3 x lO^M,^ if the particles are 
embedded in a protostellar disk of gas. Note that the inferred values of Mf^^^ and Ml°l^ are 
comparable to Jupiter and Saturn's masses and the Oort's cloud (Duncan, Quinn, & Tremaine 
1987). This mass is larger than the present mass estimate (~ O.IM^) for the Kuiper Belt today 
(Luu et al. 1997) but may be comparable to it in the past (Stern et al. 1997, Kenyon & Luu 
1999). The associated surface density of heavy elements is also comparable to that extrapolated 
from the minimum mass nebula for Rq = 70AU. 



3.4. Preserved particle scenario 

We now consider the possibility that a range of particles were formed inside a gaseous 
protostellar disk. When the surface density of the disk is depleted, particles in the outer regions of 
the disk become exposed to the stellar radiation, all particles with ad < 02 would be blown away 
by the radiation pressure. If these small particles are not the dominant IR emitters, they would 
not not need to be continually replenished by some parent particle. In their model analyses of a 
large wavelength range of IR emission, Augereau et al. (1999) adopt a distribution of particles 
with sufficiently large size to survive radiation blow outs. We follow their approach by examining 
the evolution of the surviving particles. 

Using our model parameters, we find a2 = 100 — 200^m (see Fig. |l] and Fig. |2|). (Note that 02 
does not directly depend on the gas surface density as long as ao > 02). For illustrative simplicity. 



we assume the IR-emitting particles have a characteristic size = ao = 600^m such that eg (|24|) 
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gives Mchar ~ 2.3M0. Since these particles survive the radiation blow out, they collide up to 
~ UcTi, ~ (Liiisk/ Li,){Ti,/T(i) ~ 10^ times within the life span of the central star. 

In order to determine the collisional velocity among the particles, we note that in the absence 
of any warp, the optical depth in the disk's radial direction to the host star is ~ TzR/H^. 
Since the surface brightness of the disk is spread out over an extended radial range, r,. < 1 and 
Hd/R > Ldisk/Li, for the IR emitting particles. The corresponding velocity dispersion of the 
particles would be (jrf > 3 x 10'^ cm s~^. If these particles are embedded in a gaseous background, 
Hd is expected to be < Hp which corresponds to an upper limit cr^ < lO^cm s^^. Within this 
velocity range, collisions may lead to either erosion (Hatzes et al. 1991) or coagulation if the 
particles are coated with frosts or traces of methanol as in comets Supulver [et al. 1997). Repeated 
collisions may eventually lead to an equilibrium in which fragmentation and coagulation are 
balanced to result in a size distribution similar to that in eq(|29|) (see review by Artymowicz 1997). 

In addition to the mass requirement, a population of surviving particles (with > ai and 
d-d > 0,2) need to have a spatial distribution similar to that observed ring structure whether they 
are parent particles or the dominant IR emitting particles. In the following subsection we show 
that particles which are marginally larger than 02 can congregate into rings if there are some 
non monotonic variations in the S distribution of the gas. However, if particles become larger 
than ao, they decouple from the gas flow. Nevertheless, this confinement mechanism continues to 
be effective provided dV become positive although the rate of concentration decreases with the 
particle size. Meter sized particles need ~ 10^ times longer to concentrate as they radially drift 
about lO^times slower then the particles. 

3.5. A set of quasi-static models 
3.5.1. Model: A 

We first consider the hydrostatic structure of both gas and dust disk. We set the initial 
surface density distribution of the gas to be proportional to as it is typical for protostellar 

disks. The presence of gas is invoked in order to establish a flow with velocity reversal which 
would spontaneously lead to the emergence of ring structure. In §2.1, we adopt the surface density 
distribution of the gas disk to be that in eq (1) with its mass within 70 it 20 AU to be 99M^ in 
order to augment for the solar composition. As a free parameter, we set d = 0.9 which corresponds 
to a smooth 10% surface density increase at 70 AU. The amplitude of this variation in S (see Fig. 
H and Fig. ||) is much smaller than that associated with gap formation. 

The disk structural parameters are calculated self-consistently for a dynamically stable, 
rotating, optical-thin disk. This equilibrium configuration is established on the dynamical 
time-scale r^. The increase in surface density of the gas at 70 AU corresponds to a minor peak 
in density, a local minimum in the friction time and stokes number, but it has no effect on 
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temperature and pressure scale height, because the disk is optical thin (Fig. |l]). The pressure as 
function of density also shows a small amplitude local maximum. The resultant local maximum in 
the radial pressure gradient strongly influences the local rotational profile such that the deviation 
from the Keplerian rotation speed is decreased by 50%. The velocity resulting from the Poynting 
Robertson effect is orders of magnitude smaller and it attains this local minimum as a consequence 
of the decrease in the friction time rj (see eq. 5). The last frame shows the minimum radius 
for particles to be gravitationally bound (02), which is given in eq([l8|) by the equilibrium of the 
radiation pressure, gravity and centrifugal forces. This plot also indicates the upper particle size 
limit for particles to be stopped at 70AU, which is about 700 — 800/i. 

For these disk parameters, our chosen particle's size is (a^ = ao = 600/xm) to receive maximal 
velocities and thus to show the upper limit of concentration. In Fig. ^, the solid line represents the 
resulting drift velocity of these particles due to combined effect of gas drag and radiation pressure 



in accordance to eq(15). In order to consider the relative importance of various contributions, 
we compare our main results with the dashed-line velocity curve for particles in a gravity- and 
rotation-free system, where radial radiation pressure and friction forces cancel each other. In this 
idealized case, the particles' radial velocity is positive everywhere because the friction time grows 
faster than the radiation pressure decreases with radius. The decline in Vr at 70 AU is due to the 
local minimum in the friction time-scale. 

We also show with the dotted line, the Vr distribution in a system where the radiation 
pressure is neglected and the orbital evolution is determined by gas drag only. In this limit, the 
radial drift is proportional to the friction time and the particles migrate inward everywhere. Due 
to the non monotonic variation in our prescribed S distribution near 70 AU, the deviation of the 
gas velocity from the Keplerian flow is at a minimum. This change in dV leads to a large increase 
in Vr- But with d = 0.9 and 5Ro = 20 AU in Model A, gas drag alone is not sufficient to change 
the sign of dV and Vr- In contrast, when the effects of radiation pressure and gas drag are fully 
taken into account simultaneously, a radial velocity inversion is established near 70 AU. 

For this model, particles with ad < 02 are blown away by the radiation pressure 
on the dynamical time-scale t^. When the surface density of the gas is reduced below 
S2 = 3L^/{8TrRcc^d\np/dlnR) ~ O.lg cm~^, oq becomes smaller than 02 (see eq. ^ and eq. ^) 
and there are no longer particles that drift radially outward as a result of the combined forces 
from gas and radiation. If the motion of the disk gas is only affected by a negative pressure 
gradient, all particles with < ao would be blown away by the stellar radiation pressure and 
all particles with ad > ao would become segregated from the motion of the gas, drift inwards 
with a radial velocity in accordance to eq([T7|) and a corresponding orbital evolution time-scale 
R/vr ~ {R/ Hp)'^{adPd/'^)Td ~ {R/ Hp)'^Td. In this case, all particles with sizes smaller than 
~ a2{Hp/ R)'^{Ti,/Td) would be depleted during the lifespan of HR4796A. But S2 is only marginally 
larger than the surface density of the parent particles inferred in §3.1. In §2.4, we have already 
shown that the particles' drag have a negligible effect on the motion of the gas. Eventually, when 
the surface density of the gas is reduce below Si = 3L^/(87rGM^c) ~ 10~^g cm~^, ao < ai. All 
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particles with aa < ai are blown away by the radiation pressure whereas larger particles would not 
undergo further orbital decay as the effect of head wind diminishes. 



3.5.2. Model: B 

We also computed a model with the same d but a smaller 5Rq = IQAU which represents a 
gaseous ring with a 2 times steeper surface density gradient (see Figs. ^ and ^. In this case, the 
gas rotation is hyper Keplerian with a sign change in dV such that gas drag alone may lead to an 
outward velocity near 70 AU. All particles with > 02 or > qq would migrate outwards in 
regions with a positive pressure gradient and inwards in regions with a negative pressure gradient. 
In this case there is no upper limit for the particle size in order to be stopped at 70 AU. Fig. ^ 
shows that the resulting drift velocities for Im sized objects is about three orders of magnitude 
smaller which nevertheless allows these bodies to concentrate in reasonable time in a ring at 70 
AU. 



3.6. Dynamical evolution 



The resulting radial drift velocity (Eq. can be used with the continuity equation to 
compute the global evolution of the dust surface density as a result of their orbital migration. 
In a cylindrical coordinate system, 

The local flux of dust particles is partially contributed by the systematic radial drift Vr (see Eq. 



15). Small particles with < oq, are well coupled to the disk gas. If the disk gas is turbulent, 
as in accretion disks in many astrophysical contexts (cf. Lin & Papaloizou 1996), these particles 
may undergo diffusion following the trace of turbulent eddies. For computational convenience, we 
describe the turbulent stirring as passive diffusion in the conventional Navier-Stokes equation such 
that 

h = ^dVr-^-^{R^d). (31) 

R or 

A similar approach was used by Morfill and Volk (1984) to study the reprocessing of 
meteorites in the solar nebula. Following Cuzzi et al. (1993), we assume the turbulent diffusivity 
D = u/Sc with the Schmitt number which describes the coupling between the particle and the gas 
is set to be Sc = (1 + St)^^"^ . For the gas turbulent viscosity, we adopt the ad hoc a prescription 
which is commonly used in accretion disks (Shakura & Sunyaev 1973) such that 

V = aCsHp. (32) 

The Stokes number for particles is modified from that for a laminar disk (see §2.2.1) by 

St = TfLOt- (33) 
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ujt is the circular frequency of the largest eddy (typically smaller than 0). Thus, from the definition 
of the effective viscosity it follows: 

= (34) 

p 

For the particle size we are considering, the turbulent Schmitt number is ~ 1 so that maximum 
mixing is expected to occur. 

The differential equation (]30| ) is solved numerically using a finite volume scheme. It is 
integrated explicitly in time with a simple Euler backward scheme. 



3.6.1. Model: A 

In Figure 1^ the initial radial distribution of 600/um dust particles (= lO^^S) is shown with 
the solid line. The evolution of the Ti^ distribution is shown with a series of dotted lines which 
represents time sequences separated by 10^ yr between each line. After the first 400 yr, the dust 
particles inside 66AU drift inward while that between 66 and 69 AU drift outward. The lack of 
replenishment causes interior to 66AU to decline. The results in Figure ^ indicate that Vr is 
a decreasing function of R (i.e. the magnitude of Vr is more negative at large R) inside 66 AU. 
Consequently, the depletion of the disk is more rapid at larger R and T,d becomes a decreasing 
function of R in this region. The outward drift at 66-69AU and inward drift at radius larger 
than 69AU cause the accumulation of dust particles between 67-73 AU and the depletion of dust 
exterior to 73 AU. After several 10^ yr, an equilibrium is reached in which the turbulent diffusion 
is balanced by the systematic drifts. In this equilibrium, S^i? between 67-73 AU becomes nearly 
two orders of magnitude larger than that of the disk interior to 67 AU and another two orders 
of magnitude larger than that of the disk beyond 73 AU. The actual width of the ring depends 
on the magnitude of the effective viscosity parameter. A much wider ring may be obtained for a 
larger viscosity a > 10~^. 

After the equilibrium is established, the persistence of the gas is not required. In §2.2.2, 
we show that as S decreases during the depletion of the disk, ao decreases and tj increase 
while dV is not significantly modified. While ad is still small compared with oq, the factor 
{Li^Ad/'inR'^mc + 2QdV) in eq(p!5|) remains unchanged. Although the magnitude of Vr may 
increase due to the weakening hydrodynamic drag effect, the locations for velocity reversal remain 
unaltered. Thus, dust particle concentration in the ring is not affected by the depletion of S. 
Although dust particles would generally migrate inwards after ao decreases below their size ad, 
the drag time-scale becomes so large that the radial locations of the particles are preserved. Thus, 
this scenario only needs the gas to be present for the initial few 10^ yrs in order to produce the 
ring. Afterwards it can be evaporated or be viscously diffused, while the dust distribution remains 
intact. 
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3.6.2. Model: B 

Similar model parameter are applied to a narrower ring ( with 6Rq = l^AU). In this model, 
we examine the concentration of Im-size large objects. Fig. ^ shows the same behavior like 
the smaller particles, but this time the process is about 1000 times slower. Nevertheless, the 
concentration process takes an order of magnitude less time than the life time of HR4796A. 



3.6.3. Model: C and D 

So far the upper limit for the observed gas around HR4796A is 1-7 (Greaves et al. 1999), 
which is 1-2 orders of magnitude less than we used for Model A and Model B. Now we present to 
models with only 10 of gas (see fig. ^) and 1 (see fig. P). In comparison to the first models 
the densities are smaller, the Stokes numbers for 600 /x particles are bigger and thus the resulting 
drift velocities become also smaller. There are no significant changes in the critical particle size 02- 

Thus the dynamical evolution takes 5 times longer for Model C (see fig. [lO| ) and about 25 



times longer for Model D (see fig. 11). Even the time-scales increase at lower gas densities, the 
general effect persists. 



4. Summary and Discussions 

We examined in this paper the dynamics of the dust ring around HR 4796A with two 
separate scenarios. In the first scenario, we consider the case where the dominant emitters 
of the IR radiation are /xm-size particles. Since these particles are being blown away by the 
radiation pressure of the central star, they must be replenished by large particles which are 
not affected by such a process. In a gas-free case, the parent particles must have a size greater 
than ~ 3.2 whereas in a gas-rich environment, another order of magnitude larger size is 
needed. In order to account for the magnitude of reprocess radiation, M^har ~ 4 — 8 x lO^^M® 
of /im-size particles are needed at any given time. During the life span of IIR4796A, a total mass 
^losT ~ Mchar^-k/Td ~ SOOM® is needed in the reservoir of parent particles. In order for the parent 
particles to be inconspicuous sources of IR radiation relative to the //m-size particles, a large 
fraction of Mi^gt must be contained in particles with ad > 10cm. However, the smallest surviving 
particles (with sizes ~ 100 — 200/im), could provide adequate supplies of fim size particles if a 
large fraction of their mutual collisions results in their total fragmentation. 

The main difficulty of identifying small particles as the dominant contributors of the IR 
radiation is the requirement of a non continuous bimodal size distribution. This problem would be 
removed if the particles which emit most of the IR radiation are sufficiently large (> 100/im) to 
survive the radiative blow out. In this case, a total mass < O.4M0 is needed throughout the life 
span of HR 4796A. We suggest the preservation of these modest-size particles (whether they are 
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parent particles or the dominant IR emitters) is due to the combined effect of radiation pressure 
and gas drag. 

We presented a detailed model for the formation of dust rings in circumstellar disks. Our 
model can account for rings with sharply defined edges if it is embedded in a gas disk. A very small 
amount of inhomogeneity is needed to induce the congregation of the particles. This inhomogeneity 
can be the result of a planet, but it does not necessarily has to be so. Photo-evaporation, an 
inhomogeneous initial distribution of mass after the infall, or inhomogeneities in the radial viscous 
transport of material can also account for such an anomaly in the surface density distribution. 
Thus we conclude that dust rings with sharp structures are not necessarily good indicator for 
embedded planets. 
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Fig. 1. — Model A: The radial distribution of surface density, density, temperature, friction time 
for 600 /i-size particles, the resulting Stoke-Number, pressure scale height, the speed of sound, 
the pressure, the radial pressure gradient, the deviation from the Kepler rotational profile, the 
Poynting-Robertson effect and the size (cm) of the smallest gravitationally bound particles. The 
physical parameters for this plot are given by the mass (M^ = 2.5Mq) and luminosity (L^ = 35Lq) 
of the central object. Additional model parameters include the location of the maximum in surface 
density at 70AU, the width of the inner edge as 5Ro = 20AU and the gas mass around the edge 
between Rq — 5Ro and Rq + 6Ro which is taken to be 100 earth masses. 

Fig. 2. — Model B: The radial distribution of surface density, density, temperature, friction time 
for 600 ^-size particles, the resulting Stoke-Number, pressure scale height, the speed of sound, 
the pressure, the radial pressure gradient, the deviation from the Kepler rotational profile, the 
Poynting-Robertson effect and the size of the smallest gravitationally bound particles. The physical 
parameters for this plot are given by the mass {M^ = 2.5Mq) and luminosity (L^ = SSLq) of 
the central object. Additional model parameters include the location of the maximum in surface 
density at 70AU, the width of the inner edge which is 5Ro = lOAU and the gas mass around the 
edge between Rq — 6Rq and Rq + 5Rq which is 100 earth masses. 

Fig. 3. — Radial velocities distribution of 600 /x-size particles for the model in Fig. |^. The solid 
line gives the effective radial drift. The dotted line accounts only for the gas-drag effect and the 
dashed line denotes only for the radiation pressure effect. 

Fig. 4. — Radial velocity distribution of 600 //-size particles for the model in Fig. The solid line 
gives the effective radial drift. 

Fig. 5. — Radial velocity distribution of 1 m-size particles for the model in Fig. |2|. The solid line 
gives the effective radial drift. 

Fig. 6. — Model A:Evolution of the 600/i-size particle surface density distribution for the model in 
Fig. |. The sohd line gives the initial distribution. The following lines are snapshots every 400 yrs. 

Fig. 7. — Model B:Evolution of the Im-size particle surface density distribution for the model in 
Fig. ^. The solid line gives the initial distribution. The following lines are snapshots every 10^ yrs. 

Fig. 8. — Model C: The radial distribution of surface density, density, temperature, friction time 
for 600 /i-size particles, the resulting Stoke-Number, pressure scale height, the speed of sound, 
the pressure, the radial pressure gradient, the deviation from the Kepler rotational profile, the 
Poynting-Robertson effect and the size of the smallest gravitationally bound particles. The physical 
parameters for this plot are given by the mass (M^ = 2.5Mq) and luminosity {L^, = 35Lq) of 
the central object. Additional model parameters include the location of the maximum in surface 
density at 70AU, the width of the inner edge which is 6Rq = lOAU and the gas mass around the 
edge between Rq — 6Rq and Rq + 6Rq which is 10 earth masses. 
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Fig. 9. — Model D: The radial distribution for surface density, density, temperature, friction time 
for 600 /i-size particles, the resulting Stoke-Number, pressure scale height, the speed of sound, 
the pressure, the radial pressure gradient, the deviation from the Kepler rotational profile, the 
Poynting- Robertson effect and the size of the smallest gravitationally bound particles. The physical 
parameters for this plot are given by the mass (M^ = 2.5Mq) and luminosity (L^ = SSLq) of 
the central object. Additional model parameters include the location of the maximum in surface 
density at 70AU, the width of the inner edge which is 5Ro = lOAU and the gas mass around the 
edge between Rq — 6Rq and Rq + 5Ro which is 1 earth mass. 

Fig. 10. — Model C: Evolution of the 600;U-size particle surface density distribution for the model 
in Fig. ^ The solid line gives the initial distribution. Other lines are separated by a time interval 
of 2000 yrs. 

Fig. 11. — Model D: Evolution of the 600^-particle surface density distribution for the model in 
Fig. ^ The solid line gives the initial distribution. The other lines are separated by time intervals 
of every 10^ yrs. 
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